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FIGURES  (Continued) 


Normalized  beam  spreads  for  the  long  term  and  short 
term  irradiance  profiles  for  30  = 1.0  and  6 = 0.  5 
as  determined  from  the  1/e  point  in  intensity 

The  improvement  in  on-axis  irradiance  due  to  tilt 
removal  as  a function  of  the  effective  strength  of 
turbulence 


Introduction 


Many  laser  systems  operating  in  the  atmosphere  can  be  degraded 
by  atmospheric  turbulence.  For  example,  target-illumination  systems, 
communication  systems,  radar  systems,  and  others  may  be  severely 
affected  by  turbulence  effects  which  include  beam  spreading  and  scintillation. 
Since  beam  spreading  and  its  associated  on-axis  irradiance  degradations 
are  of  practical  interest,  a detailed  theoretical  description  of  the  im- 
proved on-target  irradiance  characteristics  resulting  from  cancellation 

of  turbulence-induced  beam  wander  effects  is  desired.  Previous  treat- 
(1-4) 

ments  of  the  atmospheric  induced  laser  beam  spreading  have  been 

concerned  mainly  with  the  determination  of  the  long-term  average  beam 

spread,  i.  e.  , the  effects  of  both  beam  breathing  and  wander  are  included. 

( 5) 

A previous  semi-quantitative  analysis  indicates  that  the 
short-term  average  beam  spread  (beam  spread  with  a wander  correction) 
rapidly  approaches  the  long-term  average  beam  spread  as  the  strength 
of  turbulence  increases.  This  is  contrary  to  the  result  of  recent  un- 
published experiments^)  and  the  theoretical  analysis  presented  below. 

Rather  than  treating  the  short-term  average  beam  spread  alone,  the 
complete  analytic  expressions  for  the  short-term  average  (i.e.,  beam 
wander  cancelled)  irradiance  profile  of  a focused  laser  beam  transmitted 
through  isotropic  homogeneous  turbulence  for  a general  aperture  dis- 
tribution (and  in  particular  for  three  specified  aperture  functions)  will 
be  derived.  Then, short-term  laser  beam  spread  may  easily  be  determined 
from  the  1/e  value  in  irradiance. 

(5  7) 

Following  previous  discussions,'  * an  optical  beam  interacts 
with  atmospheric  refractive  inhomogeneities  of  all  scales.  Turbulent 
eddies  having  characteristic  dimensions  greater  than  the  beam  diameter 
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lead  to  refractive  effects  (beam  wander);  these  large-scale  inhomogeneities 
lead  to  net  wavefront  tilt.  In  this  case,  there  is  no  beam  breathing  and 
the  only  effect  is  that  the  entire  beam  dances  around  its  unperturbed 
position.  Conversely,  eddies  having  characteristic  dimensions  less  than 
the  beam  diameter  lead  to  diffractive  effects  (breathing)  about  the  instan- 
taneous center  of  energy.  The  term  long-term-average  is  used  to  imply 
averages  over  a time  long  compared  with  all  fluctuations  of  interest 
(i.  e.,  over  all  possible  ensembles  that  the  index-of-refraction  fluctuations 
can  possibly  assume).  This  means  that  the  effects  of  beam  tilt  (wander) 
and  beam  breathing  are  included  in  the  analysis  of  the  propagation  pro- 
blem when  the  long-term  average  is  taken.  Here,  however,  we  are  con- 
cerned with  the  short-term  average  irradiance  profile  about  the  instan- 
taneous center  of  energy  of  the  focused  beam.  This  case  is  of  physical 
interest  when  the  beam  dances  as  a whole  and  does  not  break  up  into 
multiple  blobs,  i.  e. , when  amplitude  fluctuations  are  small  as  in  the 
near  field  of  the  effective  coherent  transmitting  aperture.  This  is 
equivalent  to  the  requirement  that  the  propagation  distance  z « k l where 
k is  the  optical  wave  number  and  i the  smaller  of  the  initial  transmitted 
beam  diameter  (D)  or  the  long-term  spherical  wave  lateral  coherence 

length'  (PQ)»  This  is  equivalent  to  the  statement  that  the  variance  in 

7 

the  log -amplitude  of  a spherical  wave  is  much  less  than  unity  (Xg)  « 1. 
(Sharp  angular  brackets  denote  the  ensemble  average  and  we  will  drop 
the  subscript  s,  indicating  spherical  wave,  in  the  remainder  of  this 
paper.  ) Finally,  since  turbulent  eddies  of  characteristic  dimensions 
greater  than  the  beam  diameter  D are  responsible  for  producing  tilt, 
the  short-term  average  will  imply  averages  performed  over  a time 
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shorter  than  D/Vn>  where  Vr  is  the  normal  wind  velocity  (i.  e.  , assuming 


a frozen  turbulence  model).  This  indicates  that  the  short-term  average 
excludes  the  effect  of  beam  wander. 

General  considerations  regarding  the  ir radiance  profile  are  dis- 
cussed in  Section  I,  where  it  is  shown  how  the  short-term  average  irradiance 
profile  depends  on  the  modified  phase  structure  function.  In  Section  II, 
this  modified,  or  short-term  average,  structure  function  and  therefore 
the  mutual  coherence  function  (MCF)  of  a spherical  wave  is  shown  to  be 
a function  of  not  only  the  difference  coordinate  (IT),  of  two  arbitrary 
points  in  the  aperture,  but  is  also  a function  of  the  corresponding  sum  co- 
ordinate  (R).  This  is  contrary  to  the  results  employed  by  Fried'  where  sta- 
tistical independence  of  phase  with  tilt  removed  and  the  random  vector  3 defining 

tilt  was  assumed.  To  proceed  further,  a Kolmogorov  spectrum  of  tur- 

> 

bulence  is  assumed  and  three  specific  forms  for  an  aperture  function  are 
employed  to  determine  the  modified  MCF.  Using  the  modified  MCF, 
expressions  for  the  short-term  average  irradiance  profiles  are  presented 
in  Section  III.  The  aperture  functions  considered  are  uniform  with 
finite  cutoff  and  obscuration,  gaussian  with  finite  cutoff  and  obscuration, 
and  infinite  gaussian  with  semi-gaus sian  obscuration.  A numerical 
example  is  provided  for  the  case  of  the  truncated  gaussian  with  obscuration 
for  various  values  of  normalized  strength  of  turbulence  (s)  varying  be- 
tween 0.  1 and  10.  0 where  s = D/pQ.  For  comparison  the  long-term 
irradiance  profiles  will  be  exhibited  for  the  same  cases.  Plots  com- 
paring the  center  of  energy  irradiance  values  and  beam  half  width  for 
various  cases  are  also  presented. 


The  present  analysis  is  restricted  to  a weakly  inhomogeneous 
medium  with  characteristic  scale  lengths  much  greater  than  the  optical 
wavelength.  In  addition,  it  is  assumed  that  the  characteristics  of  the 
medium  do  not  change  appreciably  during  an  oscillation  period  of  the 
optical  field  at  infrared  and  visible  wavelengths.  This  condition  is  well 
satisfied  in  the  atmosphere.  The  electromagnetic  field  considered  has 
a time  dependence  given  by  the  multiplicative  factor  exp(-iirt).  The 
time-dependent  wave  equation  is  replaced  by  the  Helmholtz  equation  for 
an  inhomogeneous  medium  having  an  electrical  conductivity  and  magnetic 
permeability  equal  to  zero  and  one,  respectively.  Only  the  case  of  the 
propagation  of  a scalar  field  in  a fluctuating  medium  is  discussed.  Pro- 
files are  calculated  at  the  focal  point  with  extension  to  non-focal  points 
being  straightforward  though  nontrivial.  For  simplicity  a spatially 
homogeneous  turbulent  medium  is  assumed. 

I.  General  Considerations 

The  short-term,  long-term,  and  Fried's  approximation  " to  be 
short-term  irradiance  profile  are  obtained  from  the  extended  Huygens- 
Fresnel  principle  of  Lutomirski  and  Yura*2*.  From  this  principle,  the 
instantaneous  field  due  to  an  arbitrary  complex  aperture  distribution 
can  be  computed  by  superimposing  wavelets  that  radiate  from  all  elements 
of  the  aperture.  The  properties  of  the  medium  appear  only  in  the  propa- 
gation properties  of  the  spherical  wavelets  thus  separating  the  geometry 
of  the  problem  (aperture  distribution)  from  the  propagation  problem. 

The  instantaneous  irradiance  at  a point  p located  in  the  plane  perpendicular 
to  the  z-axis  at  a propagation  distance  z from  a transmitting  aperture  is. 

The  use  of  the  words  "Fried's  approximation"  is  used  in  the  context  that 
one  employs  the  approximations  made  in  Ref.  (7)  for  the  modulation  transfer 
function  (MTF)  to  the  MCF  used  here  to  find  the  irradiance  profile.  Fried 
never  made  beam  calculations,  he  calculated  the  MTF  in  short  exposure 
imagery. 
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w*r* 


(2 ) 

in  the  paraxial  approximation'  , 

«5>  ■ - (t*=)2JJ ^[-ik?Tl'72>]  «*  [4  ( -f  - 4)  -4 

(1) 

X exp  ^(Tj)  + **(r^)J  uA^i>  ua***2*  d2  r*^  d2?^ 

where  (r*"j ) is  the  specified  initial  complex  aperture  amplitude  distribution, 
\j(r*j ) and  i!r(r^)  are  the  perturbations  in  the  field  at  p"  due  to  the  atmosphere 
for  unit  spherical  waves  propagated  from  7^  and  F,,  f is  the  focal  range,  k is 
the  propagation  constant  2rr/\  with  \ the  wave  length,  and  for  convenience  the 
constant  E is  chosen  such  that  the  total  power  transmitted  through  the  aperture 
is  a constant  irrespective  of  obscuration  or  type  of  aperture  function. 

In  Eq.  (1)  we  have  neglected  derivatives  of  the  perturbations  iif  with 
respect  to  the  propagation  distance  as  they  are  small  quantities  compared 
to  k.^2'  After  an  average  is  performed  over  a time  that  is  long  compared 
with  all  fluctuation  periods  of  interest  (over  all  possible  ensembles  that 
the  index-of-refraction  fluctuations  can  possible  assume),  the  long-term 
(LT)  mean  irradiance  profile  [for  the  focused  case  (f  = z)]  is  obtained 


from  Eq.  (1)  as 


<I(p^)LT  = E 


ta‘ /( - [■  ^21] 


MLT(ri’r2'z) 


X uA(fpuA*(r2)  d2rx  d2r2 

where  is  the  long-term  MCF  of  a spherical  wave  located  at  the 

— — (2) 

point  p and  observed  in  the  aperture  plane  at  r^  and  r2.  ' For  the  case 

of  homogeneous  isotropic  turbulence 
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ii) 


mm 


MLT(rl’  r2'  z)  " MlT(p'z)  = exP  [-  J D^(p,z)J 
where  p is  the  magnitude  of  p*  = ^ and  D.(p,z)  is  the  wave  structure 

function  defined  as  the  sum  of  log -amplitude  and  phase  structure  functions:'81 

D*  = °X  (P'  Z)  + D0  (p'z)  (4) 

For  the  Kolmogorov  spectrum  in  the  inertial  subrange,  it  has  been  shown' 5) 
that 

D«  = 2 (p/ P0)5/3  (5) 

where  pQ  is  the  L.T  lateral  coherence  length  defined  in  Eq.  (14)  of  Ref.  1. 

(9 ) 

It  can  be  shown  that  the  phase-structure  function  DJp)  (we  drop  the 
explicit  z dependence)  can  be  written  as 

D<s=y  (p)  (p) 

where  or(p)  varies  smoothly  between  1 for  z « kp2  to  1/2  for  z » kp2  . 

In  this  paper  we  have  already  assumed  that  z « k^2  which  implies  that 
for  all  significant  values  of  p,  y(p)  = 1 and  that  the  log-amplitude 
variance  <x^ ) « 1. 

What  is  of  more  interest  here,  however,  is  the  irradiance  dis- 
tribution obtained  on  a short-term  basis  for  which  the  tilt  of  the  wave- 
front  has  been  removed,  for  example  by  reciprocity  tracking.  This  is 
accomplished  by  formulation  of  a precise  statistical  definition  of  the 
deformed  wavefront  "shape"  in  terms  of  a series  of  polynomials,  as 
defined  by  F ried.  '10*  The  instantaneous  wavefront  tilt  can  be  represented 
in  terms  of  linear  polynomials  and  this  tilt  is  removed  from  the 


wavefront  phase  by  defining 

Y(r)  = 0 (r)  - T*T 

where  Tis  a random  vector  related  to  <£(7)  in  such  a manner  that  8*7 gives 
the  best  fit  to  0(r)  in  terms  of  a least-squares  difference  over  the  aperture, 
i.  e. 

J*  d2  T W (7)  [0  (?)  - ■r’  • 7]2  = 0 

This  implies  that  the  i-th  component  of  the  vector  defining  tilt  ( i = x,  y) 
is 

fr.  <t(~r)  W (r)  d2"? 

8.  = jL  L — 1, 

1 jf  r2  w (T)  d2  r 

where  W(r)  is  a real  aperture  function  defined  such  that  W is  porportional 
to  |U^(r)  | . It  is  not  necessary  to  normalize  W as  the  normalization 

factor  cancels  in  the  definition  of  8.  (Eq.  8).  The  above  choice  for  W was 
made  to  maximize  tilt  removal.  Other  choices  were  considered  during 
this  study  but  were  found  to  give  non-optimum  results.  With  tilt  removed 
the  resulting  beam  profile  is  determined  by  the  geometry  of  the  problem 
(the  arbitrary  complex  aperture  distribution)  and  the  propagation  problem 
(diffractive  scattering  by  eddies  having  characteristic  dimensions  less 
than  the  order  of  the  beam  diameter).  The  use  of  the  extended  Huygens - 
Fresnel  principle  separates  these  two  problems. 

Averaging  Eq.  (1)  for  times  short  compared  to  D/Vn»  the  short- 
term (ST)  irradiance  profile  is  obtained  for  f = z as 
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<I(?)  >c~  = E 


(rsr)2Jj axp  [-  ~ ,(«l — ~]mst'V 


xUA,ri»  UA*(r2,d2?i  d2f2 


where 


= <exP^(ri>  + **(r2)  ‘ i®’*  (ri  ‘ r2)]^ 


is  the  "effective  short-term"  average  MCF  of  a spherical  wave  located  at 
point "p" observed  in  the  plane  of  the  transmitting  aperture. 

II.  Short  Term  Average  MCF 

To  be  general,  Eq.  (10)  is  rewritten  to  exhibit  the  dependence  on 
phase  and  log -amplitude: 

MST(r*1,  r*2)  = < exp  j x(r\)  + X(r“2)  + i[0(xr1 ) - ’ 0 * (^\  " r2  ^ I ^ (11) 

As  in  Ref.  7 it  is  assumed  that  9,  like  X ai*d  0,  is  normally  distributed. 
Furthermore,  the  distribution  of  [0(7^)  - 0fr^)  - 9 • (r£  - r2)  independent  of 
the  distribution  of  [>(1^)  + x(?^)]»  since  we  have  assumed  that  the  tur- 
bulence is  locally  isotropic  and  homogeneous.  In  contrast  to  Ref.  7,  it 
is  not  assumed  here  that  the  distribution  of  T is  independent  of  the  dis- 
tribution of  [0(7)  -T-7].  Using  conservation  of  energy,  it  is  then  easy 


to  show  that 


^T,  r2)  = exp  J - j [Dx  (p)  + D^Tj,  r2)]| 


where  D^fr^  ,r^ ) the  modified  phase  structure,  given  by 

Dv(rr1.72)=  < [ v ( r*i  ) - Y(r2)]2>  > 

= ([tffrj)  - 0<r2)]2)  + <(r*"p)2>  " 2?  * (TlXrp  - 0{TZ)]) 


I 


Substituting  Eq.  (8)  into  Eq.  (13)  and  assuming  that  the  aperture 
function  W is  circularly  symmetric,  the  modified  phase  structure  function  is 
given  by 

°v(rrr2)  = d®(p) 

2 

- (?r*7/)  W(r>W(r')  D^fT-T'l)  d2rd2r/  (14) 

+ TW(r)  (jT-  Tj)  d2r  -J r*W(r)  D^lr*-  r^|)  d2rj  , 


where 


A = ^TT^r3  W(r)  dr^ 


Equation  14  may  be  rewritten  through  the  use  of  symmetry  in  the 


angular  variable  as 


Dv(^R)  = D0(p)  P2jj(r‘  W(r)W(r')  D ,(|r-fj  ) dVd2? 

(fl  f - — _ _ 2 

+ A j-  I ri  • r W(r)  ( | r - rl  | ) d r 

r J 


{T2  ' 


l—T~  f r2  • r W(r)  D (l“- r2|)  d2r 

r. 


— ri  + r2  _ _ _ 

where  R = g and  as  above  p=  r^  - r Note  that  the  factor  multi- 

2 , 
plying  p in  Eq.  (16)  is  a constant  while  the  factors  multiplying  fT*  r^  and 

p*-*  "E,  are  functions  of  only  r^  and  7^  respectively.  The  "effective"  short- 
term structure  function  is  given  by  + A,  where  is  the  long- 

term wave  structure  function  and  A represents  the  effects  of  tilt  removal 
and  is  given  explicitly  by  the  integral  terms  in  Eq.  (16).  The  details 
of  obtaining  the  short-term  spherical  wave  MCF  for  specific  choices  of 
W is  left  to  the  Appendix. 
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III.  Irradiance  Profiles  for  Specified 
Aperture  Distributions 

In  this  section  the  expressions  for  the  short-term  irradiance  pro- 
files for  a focused  beam  transmitted  through  isotropic  homogeneous  tur- 
bulence is  presented  for  the  case  of  a truncated-obscured  uniform  amplitude 
aperture  distribution,  the  case  of  a truncated-obscured  gaussian  amplitude 
aperture  distribution  with  a gaussian  intensity  halfwidth  of  or  , and  the  case 
of  an  infinite  gaussian  with  a gaussian  halfwidth  given  by  b and  with  semi- 
gaussian  obscuration.  (The  expression  for  the  amplitude  aperture  dis- 
tribution for  this  case  is  given  by 

where  is  a constant  determined  by  the  normalization  of  the  irradiance 
profile).  The  three  cases  above  are  designated  by  the  symbol  V (=  1,2,  3). 
For  cases  v = 1,2  the  truncation  and  obscuration  radii  are  given  by  b and  a, 
respectively.  These  amplitude  aperture  functions  and  the  corresponding 
effective  short-term  MCF  are  substituted  into  Eq.  (9)  to  obtain  the 
irradiance  profiles.  To  simplify  the  expression,  a change  of  integration 
variables  in  Eq.  (9)  is  made  fromT^  and  7^  to  the  sum  and  difference  co- 
ordinates R and  "(T.  Further,  the  integrals  are  made  dimensionless  and 
dimensional  factors  are  included  into  the  constant  E.  To  accomplish  this, 

R is  divided  by  b and  "p- by  2b.  Although  the  same  symbols  will  be  used 
in  the  integration  to  determine  irradiance,  they  are  now  dimensionless. 

The  short-term  MCF  corresponding  to  the  three  aperture  distributions 
is  given  in  terms  of  the  dimensionless  p”,  £ from  Eqs.  (12)  and  (16)  by 
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mst v<P’R>  = exp I ' e (p5^3  + Cv  p2) 

,e\  5/3 

+ (-f)  17-  ri  Gv(rl>  - P‘  ^ GvM| 

where  e = 2b/pQ,  Tj  = R + p-  and  = R - p\  Cv  and  Gv(x)  are  given  in 
the  Appendix,  and  the  Kolmogorov  spectrum  of  turbulence  is  assumed. 

Note  that  for  the  Kolmogorov  spectrum,  the  long-term  average  MCF  is 
obtained  from  Eq.  (17)  by  setting  both  C v and  0 while  Fried's  approxi 

mation^  ^ to  the  short-term  average  MCF  may  be  obtained  from  Eq.  (17) 
by  changing  the  sign  of  Cv  and  setting  Gv  = 0. 

The  expressions  for  the  short-term  average  focused  ir radiance 
profiles  for  cases  V = 1,2  are  obtained  from  Eqs.  (9)  and  (17)  in  a manner 
similar  to  that  used  in  the  Appendix  for  obtaining  the  quantity  and  are 
presented  without  further  comment 

<iv(p))st  = — (-inr)  |^-  sv 

Hv  (p-,  R - ~p)  RdR 
Hv  (-fT,  It  - -p)  RdR 
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9 = a-9,  6 = a/b  , 


3q  = b/o<o  is  the  ratio  of  the  aperture  radius  to  the  transmitted  gaussian 
intensity  halfwidth,  P is  the  total  power  transmitted  through  the  aperture, 
Gp  Cj,  G2>  and  C2  are  defined  in  Eqs.  (A. 8),  (A.  17),  (A. 20),  (A. 21), 
respectively  and  the  angle  between  ~p'  and  R is  denoted  by  9. 

For  the  case  of  the  infinite  gaussian  (V  = 3),  the  short-term 
irradiance  profile  is  given  by 


and  Cj,  Gj,c  and  e are  defined  in  Eqs,  (A.  29),  (A. 30),  and  (A. 31) 

respectively.  7 

2 c 

Expressions  (18)  and  (19)  when  normalized  by  (-r ) 

kb  nb 

and  evaluated  as  a function  of  — are  only  dependent  on  three  parameters. 

z 

These  are  the  strength  of  turbulence  s,  the  obscuration  ratio  6,  and  the 
taper  ratio  9Q.  For  the  non-focus ed  case  this  would  no  longer  be  true. 

TT 

Indeed,  the  results  would  also  depend  on  the  Fresnel  number  — r and 

K Z 

the  ratio  of  the  focal  range  to  the  propagation  range  — . 

z 


IV.  Numerical  Results 


The  expression  for  the  short-term  focused  irradiance  profile 

, , 2. 


normalized  by 


(#) 


for  cases  v = 1,  2 have  been  evaluated 


numerically.  A representative  example  of  the  short-term  irradiance  pro- 


file for  the  truncated  obscured  gaussian  aperture  distribution  (V  = 2)  is 

presented  in  Fig.  1 for  6 = a/b  = 0.  5 and  9 = h/a  = 1.  The  curves  in 

00 

this  figure  are  shown  for  e = 2b/  pQ  = 0.  1,  0.  25,  0.  5,  1.0,  2.  5,  5.  0, 

7.5,  and  10.  0.  For  comparison  we  have  plotted,  in  Fig.  2,  the  long- 
term focused  irradiance  profiles  for  the  same  values  of  parameters  used 
in  Fig.  1.  From  these  and  other  calculated  profiles  the  on-axis  irradiance 
and  beam  spread  may  be  determined  as  a function  of  the  strength  of 
turbulence.  The  on-axis  irradiance  values  computed  from  the  short-term, 
long-term,  and  Fried's  short-term  expressions  are  shown  in  Figs.  3,  4, 
and  5 as  a function  of  s.  [in  Fig.  3,  a uniform  amplitude  distribution 
(V  = 1)  with  an  obscuration  of  0.  5 was  used  to  calculate  the  three  on-axis 
irradiance  curves  shown.  A truncated  gaussian  amplitude  distribution 


(v  = 2)  with  an  obscuration  of  0.  5 and  a taper  ratio  of  6q  = 1 was  used 
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Normalized  long  term  irradiance  profiles  for  various 
strength  of  turbulence  values. 


in  calculating  the  on-axis  irradiances  shown  in  Fig.  4 while  for  Fig.  5 
a truncated  gaussian  with  no  obscuration  and  a gaussian  intensity  half- 
width equal  to  one -fourth  the  aperture  radius  (i.  e.  , = 4)  was  used]. 

This  last  case  is  similar  to  what  would  have  been  obtained  for  an  infinite 
gaussian  distribution  with  no  obscuration.  Note  that  Fig.  4 corresponds 
to  the  cases  presented  in  Figs  1 and  2.  In  Fig.  6 the  short-term  and 
long-term  turbulence  induced  beam  spreads  normalized  to  the  s = 0 
values  (i.  e.  , vacuum)  are  plotted  as  a function  of  s for  the  same  cases 
considered  in  Figs.  1 and  2.  The  beam  spreads  are  determined  from 
1/e  of  the  maximum  value  of  the  irradiance  curves.  Finally,  in  Fig.  7 
we  plot  the  improvement  in  on-axis  irradiance  due  to  tilt  removal 
(IsT^LT^  aS  a ^uncti°n  the  effective  strength  of  turbulence  s ^ for 
several  different  cases  of  6 and  9q.  The  effective  strength  of  turbulence 

is  defined  as  e for  V = 1 while  for  V = 2 e r , = D /p  for  9 > 1 where 

eff  eff  o o 

the  effective  aperture  diameter  is  that  given  by  Kerr  ^ 1 1 ' 


From  an  analysis  of  the  calculational  results  (see  Fig.  3 - 5), 
it  was  found  that  Fried's  approximation  to  the  short-term  MCF  may  be 
used  to  determine  the  irradiance  profiles  for  the  case  of  a uniform 
aperture  distribution  (V  = 1).  Excellent  agreement  was  obtained  when 
using  Fried's  MCF  compared  to  that  using  the  more  exact  short-term 
MCF  for  all  cases  considered,  i.  e.  6 £ 0.  5,  g < 10.  On  the  other 
hand,  it  was  found  that  the  use  of  Fried's  MCF  led  to  significant  errors 
in  the  irradiance  profiles  for  the  truncated  gaussian  aperture  distribution 
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The  improvement  in  on-axis  ir radiance  due  to  tilt  removal 
as  a function  of  the  effective  strength  of  turbulence. 


(V  = 2)  when  c > 2 and  6q>  2.  In  fact  for  3q>  2 and  e > 7,  the  use  of 
Fried's  MCF  leads  to  negative  irradiance  values. 

In  Fig.  6 it  is  seen  that  long-term  beam  spread  increases  in  a 
linear  manner  with  s with  a break  occurring  near  e = 7.  This  change  in 
slope  is  due  to  the  large  obscuration  used  for  the  case  presented.  The  linear 
dependence  of  long-term  beam  spread  has  been  noted  previously^. 

On  the  other  hand,  the  short-term  average  beam  spread  remains  con- 
stant, then  increases  slowly  for  e Z 5.  In  fact,  if  a different  definition 
of  beam  spread  were  used  (i.  e. , the  position  of  the  first  minimum)  the 
short-term  beam  spread  would  remain  constant.  This  short-term 
beam  spread  dependence  on  e has  been  seen  for  all  cases  considered 
and  is  in  contrast  to  Ref.  5 where  a minimum  in  the  ratio  of  short  to 
long-term  beam  spreads  is  exhibited.  In  Ref.  5,  the  short-term  beam 
spread  was  determined  from  the  short-term  lateral  coherence  length 
obtained  from  Fried's  approximation  to  the  short-term  MCF.  In  con- 
trast to  this  previous  work,  we  find  that  it  is  in  general  not  possible 
to  characterize  the  short-term  beam  spread  by  an  effective  short-term 
lateral  coherence  length  (see  below).  We  also  find  the  recent  phenomeno- 
logical analysis  of  beam  spread  presented  by  Kerr^^  to  be  somewhat  in 
contrast  with  the  results  obtained  here.  Indeed,  Kerr's  results  imply 
larger  on-axis  irradiance  values  and  smaller  beam  spread  angles  for  e < 1.69 
than  for  the  turbulence  free  case.  ' Assuming  that  Kerr's  results  apply 
only  for  c j£1.7,  we  are  lead  to  beam  spreads  larger  (e.  g.  , factors  of  two 
for  large  values  of  e)  than  found  in  this  paper.  The  present  analysis  of 
beam  spread  indicates  that  the  beam  center  retains  its  shape  and  that 
the  excess  energy  has  been  scattered  into  large  angles  by  the  smaller 


. ..  --- 
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scales  of  turbulence.  The  results  then  look  like  the  turbulence  free 

irradiance  decreased  by  an  extinction  coefficient  plus  the  scattered 

irradiance  which  has  only  a weak  angular  dependence.  Therefore  as  long 

as  the  peak  extincted  irradiance  is  significantly  larger  than  the  peak 

scattered  irradiance,  it  is  not  possible  to  express  the  short-term  beam 

ST  ST 

spread  in  terms  of  z/kpQ  , where  pQ  is  a so-called  short-term 

(5) 

lateral  coherence  length. 

From  Fig.  7 it  is  seen  that  the  irradiance  improvement  obtained 
by  tilt  removal  is  maximized  for  7.  5 < ~e££  - 10  for  the  cases  considered.^ 
The  maximum  improvement  found  was  approximately  5.  62  for  V = 1. 

This  improvement  decreased  slowly  as  increased.  Improvement  in 
irradiance  for  the  obscured  aperture  is  less  than  that  for  the  unobscured 
cases  with  maximum  improvement  occurring  at  smaller  values.  This 

is  not  surprising  for  large  obscurations  since  a large  portion  of  the  co- 
herent aperture  has  been  degraded. 

In  this  paper  we  have  written  down  an  expression  for  the  focused 

short-term  irradiance  profile  for  a beam  transmitted  through  a turbulent 

2 

atmosphere  which  is  valid  for  a propagation  range  z < kf.  . (Recall  that  l 

is  the  smaller  of  the  initial  beam  diameter  D or  the  long-term  spherical 

wave  coherence  length.  ) This  expression  is  more  exact  but  also  it  is  more 

difficult  to  calculate  numerically  than  the  results  obtained  by  using  Fried's 
(51 

approximation.  ' ' Fortunately  it  has  been  found  that  P'ried's  approximation 
can  be  used  with  good  results  to  determine  the  MCF  and  therefore  irradiance 
profiles  for  weak  to  moderate  turbulence,  e < 5,  and  for  gaussian  aperture 
distributions  that  are  strongly  truncated,  3q~2. 

"tNote  that  B = 0 indicates  the  uniform- illumination  case.  The  existence  of 
a maximum  in  Fig.  7 does  not  imply  the  existence  of  a minimum  in  the  ratio 
of  short  term  to  long  term  beam  spreads  using  the  definition  of  beam  spreads 
given  above. 
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1 

Finally,  we  have  noted  that  the  short-term  average  beam  spread 

does  not  increase  with  e but  remains  constant.  Obviously  this  does  not 

2 

hold  as  s increases  without  bound,  since  we  have  specified  that  z < kf  . 

Setting  z = kpQ  , where  pQ  = (.  545  Cn^k^z)  ^ ^ (Ref.  5),  implies  that 
• 545  C = (k^^z**^)  or  p = (z/k)*^,  where  C ^ is  the  index  of 
refraction  structure  constant.  This  means  that  if  the  turbulence  becomes 
so  strong  (large  cj^)  that  pQ  becomes  smaller  than  the  first  Fresnel  zone 
for  the  propagation  range  z or  s'  > kD^/z  (the  Fresnel  number)  than  the 
above  theory  does  not  hold.  When  s > kD  /z  theno’(p)  (see  Eq.  5 and 
following  discussion)  can  no  longer  be  set  to  unity  and  must  be  included 
rigorously  in  the  determination  of  D^. 


APPENDIX  A 


Calculation  of  the  Modified  Phase  Structure 
Function  for  Three  Specific  Aperture  Functions 


In  this  appendix,  Eq.  (1)  is  obtained  from  Eqs.  (12)  and  (16)  for 
the  Kolmogorov  spectrum  and  three  amplitude  aperture  functions  dis- 
cussed in  the  text. 

Since  we  have  assumed  that  z « kp  , the  phase  structure  function 
becomes 

D?(P)-2(p/pol5/3  (A.  1) 

where  pQ  is  the  long-term  lateral  coherence  of  a spherical  wave  length^ \ 
When  p is  made  dimensionless  as  discussed  Section  III,  Eq.  A.  1 becomes 


(A.  2) 


n v _ ->  5/3  5/3 

D0(p)  - 2 e p 

We  consider  three  specific  cases. 

1)  Uniform  Aperture  Amplitude  Distribution  With  Truncation 
And  Obscuration 

Here 


W^r)  =1  a < | r | < b 

= 0 0 £ |r  | < a (A.  3) 

I r | > b 

where  "b"  is  the  aperture  (truncation)  radius,  and  a the  obscuration  radius. 
First,  Eq.  (16)  is  rewritten  as 


d (p;  r)  = D^(P)  - Xj(p)  + ya  nr.  Fj)  - y1  (p  r2)  , 


(A.  4) 
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) is  the  second  and  third 


where  "r^  ^ = & i ~2~  an<^  ( P)»  Y j(  p , r ^ £ 
term  of  Eq.  (16),  respectively.  Substituting  Eq.  (A.  3) into  Eq.  (15) 

yields 

. _ __4_ 

nb‘ 


L1 


-4  (1  - 64) 


-1 


(A.  5) 


where 


6 = b/a 


Upon  substituting  Eq.  (A.  1)  and  (A.  3)  into  (A.  4)  becomes 


(A.  6) 


P o rl  0 a 

Upon  performing  the  angular  integration,  using  the  series  expansion  for 
the  Hyper  geometric  function^  12\  and  expressing  the  results  in  terms  of 
dimensionless  variables,  Eq.  (A.  6)  becomes 

5/3 


Yj(T.  *i>  = -2ff)  p*  7i  Gi<rr6> 


(A.  7) 


where 


c>,rl'  6)=  T (‘-*4)'1/  4l'31777V)dr  ,A-8' 

and  F is  the  Hyper  geometric  function  of  the  indicated  variables.  To  proceed 
further,  Eq.  (A.8)mustbe  evaluated  nume ric ally 

From  Eqs.  (16),  (A.  1),  (A.  3),  and  (A.  4)  we  obtain 


X,(p)  = 


. 2 2 2tt  b 2tt  b 

A1  P 


173 


////  T • r/|'r’-r/|^^  rdr  r7  dr7  d©  d©7  (A.  9) 


0 a 0 a 
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To  evaluate  this  fourth  degree  integral  a change  of  variables  is 
made  from  7,  7'  to  the  sum  and  difference  coordinates  Tf  = r ^ r and 
5 = r - ~r'  respectively.  The  range  of  integration  of  § is  over  a circle 
of  radius  2b  while  the  range  of  integration  of  If  is  over  the  overlap  area 
of  two  circles  of  radius  "b"  and  obscuration  radius  "a"  whose  centers  are 
separated  by  a distance  of  §.  Since  the  amount  of  overlap  and  obscuration 
changes  with  §,  Eq.  (A.  9)  is  broken  into  four  terms  to  account  for  the 
differing  dependence  on  §.  These  represent:  1)  For  0 < % < 2b,  the 
overlap  area  of  two  circles  of  radius  b separated  by  a distance  2)  For 
0 < | < b + a,  the  removal  of  the  obscuration  areas  from  the  two  circles 
of  radius  b (for  § >b  + a there  is  no  obscuration  within  the  overlap  area 
of  two  circles);  3)  For  0 < § < 2a,  the  addition  of  the  overlap  area  of  the 
obscuration  circles  of  radius*a"since  this  area  was  removed  twice 
(incorrectly)  in  2 (there  is  no  overlap  of  the  obscuration  areas  for  5 > 2a; 
and  4)  For  b-a<?<b+a,  the  addition  of  the  obscuration  areas  outside 
of  the  overlap  area  of  the  circles  of  radius  b.  Figure  A-  1 is  a geometric 
representation  of  the  integration  area  for  the  T variable.  Note  that  in 
Fig.  A.  1-b  the  obscuration  area  is  shown  non -over lapped  for  clarity. 


w 


From  the  above  discussion  and  upon  performing  the  angular  integration  in 
§ immediately,  Xj  can  be  expressed  in  terms  of  dimensionless  variables 
(§  and  T)  divided  by  2b  and  b,  respectively)  as 


Xj  = 64  nb  Aj  (e) 
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(I,  + I_  + I,  +lA) 


(A.  10) 


where  the  L correspond  to  the  terms  discussed  above  and  the  corresponding 
geometry  is  exhibited  in  Fig.  A-l.  The  integration  is  over  the  double 
crossed  hatched  area  in  Fig.  A-l  through  symmetry  of  the  angular  variables. 
The  relevant  integrals  L are  given  by 

-1-  1 


cos 


[1  = J 58/3d5  J d0  J 


(T  - 2 Mcos0)  Tld T) 


§ sec0 


iMt) 


(A.  11) 


180 


47311 


where  by  a change  in  variables  T|  is  now  centered  in  the  left-hand  circle  of 
Fig.  A-l  a and  T represents  the  gamma  function.  The  second  integral  appear- 
ing in  Eq.  (A.  10)  is  given  by 
1+6 


I. 


11  O 

j S8/3dsj*  d0  f 


T)dT)  (T)  + 2PTlcos0) 


3tt  4/1  + 6\ 
44  6 \ 2 / 


0 0 
11/3 


(A.  12) 


where  the  variable  T)  is  centered  in  the  right-hand  circle  of  Fig.  A-l.b 
and  the  integration  is  performed  over  the  top  right  half  obscuration  area. 
The  third  integral  in  (A.  10)  is  similar  to  the  first  except  that  0<  § £ 6.  This 
integral  is 


da*. 
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given  by 


6 cos" * |/6  6 

1 3 =J?8/3  <**  f d9J  (Tl2  - 2Tl?co«e)TldTl 

0 0 ?cos0 

(A.  13) 

=623/3I1 


In  order  to  analytically  reduce  the  integral  occurring  in  1^,  the 
integration  area  in  T)  (Fig.  A-l  d)  is  broken  into  two  terms  described 
schematically  in  Fig.  A-2. 

(a)  AREA  TO  BE  (b)  AREA  IN  lc|  AREA  BETWEEN 

CONSIDERED  SEGMENT  EDGE  OF  APERTURE 

OF  CIRCLE  AND  SEGMENT 

BOUNDARIES 


Fig.  A-2  Schematic  of  the  Integration  Area  for  the  T] 
Variable  in  Fig.  A- Id 


Using  this  division  1^  is  written  as 
1+6 
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Of  = cos 
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The  first  term  in  the  integration  over  T]  is  shown  schematically 
in  Fig.  A-2b  (the  origin  of  the  11  integration  is  at  the  center  of  the 
obscuration  circle)  and  the  second  term  is  represented  by  Fig.  A-2c, 
(the  origin  of  the  T)  integration  is  at  the  center  of  the  large  circle). 

On  performing  the  inner  integrations,  Eq.  (A.  14)  becomes 
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and  9 =o-0.  To  proceed  further,  Eq.  (A.  15)  must  be  computed  numeri- 
cally. In  terms  of  the  dimensionless  variable  p,  X(p)  is  then  written  as 
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X( p)  = -pCj  s p 
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where 
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From  Eqs.  (5),  (12),  (16),  (A.  1),  (A.  3),  (A.  4),  (A.  7),  and  (A.  17)  the  short- 
term MCF  for  the  uniform  amplitude  aperture  distribution,  in  terms  of 
dimensionless  variables,  is  given  by 


MgT  = exp 


5/3  . 5/3  , n 2. 

: ( p + -Cj  p ) 


5/3 


+ fr)  tP*  G^rj)  - p • r2  Gj  (r2)] 
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(7) 

When  Fried's  approximation  is  made  , the  sign  of  is  changed 
and  Gj  is  set  to  zero.  However,  the  resulting  expression  does  not  agree 
with  the  results  found  in  Eq.  (10)  Ref.  5.  Reference  5 appears  to  be  in- 
correct due  to  two  errors.  The  first  is  a calculational  error  in  determining 
Eq.  (5)  Ref.  5 and  the  second  is  the  replacement  of  p by  pt  and  integration 
of  t from  0 to  1 in  the  determination  of  tilt  removal  to  account  for  the 
difference  between  plane  and  spherical  wave  propagation.  This  effect  had 
already  be  accounted  for  in  the  long-term  phase  structure  function  and  was 
included  numerically  into  the  long  term  phase  coherence  length  pQ.  This 
can  also  be  seen  by  referring  to  Ref.  4,  Eqs.  (26)  - (28). 

2)  Gaussian  Aperture  Amplitude  Distribution  with  Truncation 
and  Obscuration 
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For  this  case  the  unnormalized  aperture  function  is  given  by 


2,  2 

= e"r  ao  a < | r | < b 

= 0 otherwise 

where  is  a measure  of  the  1/e  point  in  irradiance.  The  normalization 
o 

" A^"  is  obtained  from  Eqs.  (1  5)  and  (A.  19)  as 
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and  the  symbol  8q  is  used  to  represent  b/aQ.  Without  going  into  as  much 
explanation  as  was  used  in  presenting  Eq.  (A. 7)  and  (A.  16)  for  the  uniform 
case,  we  write 
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The  L terms  represent  integrals  for  which  the  integration  area  for 
the  T)  variable  is  the  same  as  that  discussed  above  for  Eqs.  (A.  11)  - (A.  15) 
and  Fig.  A.  1).  Different  choices  are  made  for  the  order  of  integration  to 
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Equation  (A. 22)  is  equivalent  to  Eq.  (A.  11)  though  an  integration  over  T 

2 ~ * 

has  been  performed  and  a factor  of  -(4  3q  ) removed  from  the  integral. 
Further,  the  origin  of  the  T)  integration  is  at  the  center  of  the  overlap  area 
of  the  two  large  circles  in  Fig.  A -la.  The  integral  has  an  1)  integration 
over  the  overlap  area  of  the  obscuration  and  is  found  in  a manner  similar 
to  that  used  in  obtaining  Eq.  (A.  22)  (See  Eq.  (A.  13)  and  Fig.  A-lc).  The 
results  are 
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where 


R2  = -l COS0  + 


f~2.  P2  . 2fl 

V 6 - s sin  u 


The  integral  is  calculated  in  a manner  similar  to  that  used  in  calculating 
Eq.  (A.  12)  except  that  here  a different  aperture  function  W has  been  used 
and  that  the  dimensionless  variable  I is  now  integrated  between  0 and 
y 1 + 36^  fl.  The  reason  for  this  is  that  it  is  difficult  to  express  an  integral 
for  the  obscuration  area  outside  of  the  overlap  area  of  the  large  circles 
which  is  reducible  analytically  for  5 greater  than  yr  + 362^  . This  upper 
limit  in  I is  determined  by  the  intersection  of  the  large  circle  with  the 
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obscuration  circle  such  that  T)  (from  the  center  of  the  overlap  area)  is 
tangent  to  the  obscuration  circle  at  that  point.  With  the  angular  integration 
over  1)  performed  (here  from  the  center  of  the  obscuration  circle)  we 
obtain 
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where  I (x)  = e I.  (x)  and  I.  (x)  is  modified  Bessel  function  of  order  i, 


(12) 


(Do  not  confuse  L (x)  with  symbols  L above).  The  integral  1^  represents 
the  same  quantity  calculated  in  Eqs.  (A.  14)  and  (A.  15)  except  that  the  upper 
limit  on  5 is  Jl  + zhz/z  and  that  the  origin  of  the  T]  integration  is  the  center 
of  the  overlap  area.  Upon  performing  the  T]  integration  we  obtain 
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where  again  a factor  of  ( 4B  ) has  been  removed  and  included  in  the  constant 


terms  of  Eq.  (A.  21)  and 
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Finally,  for  >/ 1 + 362 /2  < § < (1  + 6)/2  the  integral  1^  corresponds  to  that 
part  of  the  obscuration  area  inside  the  overlap  of  the  large  circles.  The 
resulting  integral,  upon  performing  the  T)  integration,  is 
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where  = -l^* 

The  MCF  for  this  case  is  obtained  from  Eq.  (A.  18)  by  replacing 
Gj,  Cj  by  G2  and  C7,  respectively. 

3)  Infinite  Gaussian  Aperture  Amplitude  Distribution  with 
Semi -Gaussian  Obscuration 


The  appropriate  unnormalized  aperture  function  that  was  used  is 


given  by 


2 ..Z  / 2,  2\2 

= e"r  /b  (l  - e"r  U / , 


where  b represent  the  gaussian  halfwidth. 

This  is  an  unusual  aperture  function  and  was  chosen  so  that  the 
expression  for  the  irradiance  profile  (Eq.  19)  is  easily  expressible  in 
terms  of  fT,  R . For  this  aperture  function  the  determination  of  the  modified 
phase  structure  function  is  simplified  since  the  integrations  occurring  in 
Eq.  (A.  4)  are  carried  out  over  all  space.  The  appropriate  normalization 
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constant  A^,  and  are  listed  below 
A3  = 2 [l  ’ 2 (f ) + (f ) 


/(TTb4) 
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where  M(a,b,x)  is  the  confluent  hype rgeometric  function'  and 
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In  the  above 
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1 1,1  1 1.1  1 1.2  1 1.3 

~t  - tj  ’ ~7 = T = rr  + ~i’  7 = ;r  + 7T  ■ 

c b adb2aebafb2a 
1 _ 1 d2  1 _ 1 f2 

g2  7 4a4  ’ 7 7 4a4 

Again,  C^,  are  substituted  for  C^,  G^  in  Eq.  (A.  18)  to  obtain  the  short 
term  MCF. 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operation*  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigation*  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  system*.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
Laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics*  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high- power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  ot  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  frMcture  mechanics  to  stress  cor- 
rosion and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  tlie  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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